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Abstract 

Higher-order corrections of Einstein-Hilbert action of general relativity can be 
recovered by imposing the existence of a Noether symmetry to a class of theories 
of gravity where Ricci scalar R and its d'Alembertian OR are present. In several 
cases, it is possible to get exact cosmological solutions or, at least, to simplify 
dynamics by recovering constants of motion. The main result is that a Noether 
vector seems to rule the presence of higher-order corrections of gravitaty. 
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1 Introduction 

The issue of recovering a suitable effective action seems to be one of the main ways to 
construct a consistent theory of quantum gravity [IJ. Starting from pioneering works 
of Sakharov ||, the effects of vacuum polarization on the gravitational constant, i.e. 
the fact that gravitational constant can be induced by vacuum polarization, have been 
extensively investigated. All these attempts led to take into account gravitational actions 
extended beyond the simple Einstein-Hilbert action of general relativity which is linear 
in the Ricci scalar R. 

At the beginning, the motivation was to investigate alternative theories in order to 
see if gravitational effects could be recovered in any case. 

The Brans-Dicke approach is one of this attempt which, asking for dynamically induc- 
ing the gravitational coupling by a scalar field, is more coherent with the Mach principle 
requests [§]. 

Besides, it has been realized that corrective terms are inescapable if we want to obtain 
the effective action of quantum gravity on scales close to the Planck length (see e.g. (4J). 
In other words, it seems that, in order to construct a renormalizable theory of gravity, we 
need higher-order terms of curvature invariants such as R 2 , R^ v R^ u , R^ ualS R^ ua/ 3, ROR, 
R\3 k R or nonminimally coupled terms between scalar fields and geometry as (p 2 R. 

Stelle |5j constructed a renormalizable theory of gravity by introducing quadratic 
terms in curvature invariants. Barth and Christensen gave a detailed analysis of the 
one-loop divergences of fourth-order gravity theories giving the first general scheme of 
quantization of higher-order theories |], [7j. Several results followed and today it is well 
known that a renormalizable theory of gravity is obtained, at one-loop level, if at least 
quadratic terms in the Riemann curvature tensor and its contractions are introduced [|TJ. 
Any action, where a finite number of terms involving power laws of curvature tensor or 
its derivatives appears, is a low-energy approximation to some fundamental theory of 
gravity which, up to now, is not available. For example, string theory or supergravity 
present low-energy effective actions where higher-order or nonminimally coupled terms 
appear. ||. 

However, if Lagrangians with higher-order terms or arbitrary derivatives in curvature 
invariants are considered, they are expected to be non-local and give rise to some char- 
acteristic length /o of the order of Planck length. The expansion in terms of R and DR, 
for example, at scales larger than Z produces infinite series which should break near / 

With these facts in mind, taking into account such Lagrangians, means to make 
further steps toward a complete renormalizable theory of gravity. For instance, Vilkovisky 
Q considered a non-local Lagrangian of the form i?/(D)i?, where 

in order to construct an effective action of quantum gravity. Immediately, one realizes 
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that it can be approximated by the sum 



1.21 



i=0 



1 



so that we get (2k + 4)-order field equations. Also the case has to be mentioned 



since it can be regarded as a conformal anomaly |jTU |. 

We have to do an important remark at this point. In the effective action, terms as 
Rin l R have to be taken into account since pure terms as O l R are total divergences and 



can be ignored. Terms as □*i?n J i? can be integrated by parts giving Rn l+3 R |TT|, O, [13 



We are restricting to Lagrangians containing Ricci scalar and its derivative since in 
this case it is quite straightforward to obtain conformal transformations relating higher- 



order gravity theories to general relativity with a certain number of scalar fields fL4 



Furthermore, only these Lagrangians are interesting for constructing cosmological 
models (at least if we restrict the discussion to homogeneous cosmological models). 

However, we have to consider the fact that in the limit of the classical theory, higher- 
order Lagrangians give rise to superfluous degrees of freedom. This is a controversy in 
literature [T^, [T7|], which, in our knowledge is not solved. Some authors [0, discuss 
the possibility that only solutions which are expandable in powers of h are self-consistent, 
others [19|] consider such superfluous degrees of freedom as phases of oscillations around 
the Friedman behaviour |[20|| . In any case, the transition to our classical observed universe 
has to be accurately discussed. 

In this paper, we are going to discuss if such extra-terms in the effective gravitational 
action can be recovered by asking for symmetries of a Lagrangian which generic form is 



-gF(R, OR). 



;i.s) 



We are using the so called Noether Symmetry Approach which was extensively used to 
study nonminimally coupled theories of the form 



C 



-fj 



F(tp)R+-V^V^ 



V{ip) 



;i-4) 



where V M is the covariant derivative. In [21], ^2], [23], [21 |, it was shown that asking for the 
existence of a Noether symmetry 



L X C = 0^XC = 



;i.5) 



where Lx is the Lie derivative with respect to the Noether vector X, it is possible to 
select physically interesting forms of the interaction potential V((p) and the gravitational 
coupling F(<p). The scalar field p is generic and it can represent the matter counterpart 
in an early universe dynamics. 

The existence of Noether symmetries allows to select constants of motion so that the 
dynamics results simplified. Often such a dynamics is exactly solvable by a straightfor- 
ward change of variables where a cyclic one is present. 
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Here we want to apply the same method to higher-order theories of the form ([L 
In particular, we take into account suitable minisuperspaces whose degrees of freedom 
are a, the scalar factor of the universe, R, the Ricci scalar, and OR the d'Alembertian of 
Ricci scalar, related among them by some Lagrange constraints. As we shall see in next 
section, using a Friedman-Robertson- Walker (FRW) metric it is possible to reduce the 
Lagrangian Q1.3Q to a point-like one and then apply the Noether technique. The main 
result is that several fourth, sixth and eight-order interesting Lagrangians are recovered 
by asking for the Noether symmetry. For example, a term as V ROR, which is a part of 



the so-called 03 anomaly pq , |26| is connected to the existence of the Noether symmetry. 

A similar result works for F R 3 ^ 2 which is related to the Liouville field theory by a 
conformal transformation. The same technique gives Lagrangians like 



C = y^(F 1 R + F 2 R 2 



or 



£ = (FiR + F 2 R 2 + F 3 ROR) , 



where are constants, widely studied in literature, e.g. [|17], [27J |28], |29], [3(J, ^TJ. In 
conclusion, it seems that the Noether approach is related to the recovering of one-loop 
and trace anomaly corrections of quantum gravity. 

The paper is organized as follows. Sect. 2 is devoted to the discussion of the generic 



point-like Lagrangian which can be recovered from theories like ( |1.3| ). In Sect. 3, we apply 
Noether vector to such a Lagrangian showing that the existence of a symmetry selects 
its form. Fourth-order and higher than fourth-order models are discussed in Sect. 4 and 
5, respectively. In both cases we study the related FRW cosmology. Discussion and 
conclusions are drawn in Sect. 6. 

2 Higher-Order Point-like Lagrangians and Equations of Motion 

A generic higher-order theory in four dimensions can be described by the action 

A = J d A x^^F{R, DR, D 2 R, . . . , D k R) . (2.1) 

We are using physical units 8ttGn = c = h = 1. Equations of motion can be deduced by 
the method worked out in [|T^, |S2| 



GT = g{-g> lu (F-gR) + (g» x g™-g^g x °)g ;Xa + 

+ ^EE(^V CT + ^V CT )(n^) ; .fn J 

A i=l 7=1 V 



dn*R 



- g" v g Xa 



-In 



dF 
dn*R 



where 



QP» = R H» _ -g^R. 
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(2.2) 



(2.3) 



3 



is the Einstein tensor, and 



h \ dU3R ) 



(2.4) 



As we said, these are pure gravity (2k + 4)-order field equations. Matter can be taken 
into account by introducing the stress-energy tensor T^ v of a (non) minimally coupled 
scalar field UTTl |33|. 



For the sake of simplicity, let us restrict to the Lagrangian ( |1.3| ). In this case, we 
have eight-order field equations which becomes of sixth-order if the theories is linear in 
□i?. To apply Noether Symmetry Approach, let us take into account the point-like FRW 
Lagrangian 

C = C(a, a, R, R, DR, (OH)) . (2.5) 

The variables R and OR can be considered independent and, by the method of Lagrange 
multipliers, we can eliminate higher than one time derivatives (for fourth-order case see 
e.g. f34|). The action related to Lagrangian (|1.3| ) becomes 



A = 2vr 2 dtl a 3 F - \ 1 



R+61- + 
\ a 



2 k 

+ ~ 

a 2 



-A, 



DR-R-3 



R 



Ai,2 are given by varying the action with respect to R and OR, that is 

. 3 dF 3 dF 

— Q> n ~ i ^2 — O tt, — • 

dR' d(nR) 
After an integration by parts, the (Helmholtz type) point-like Lagrangian is 



,dF_ 
dR~ 



C = 6aa 2 ^—+6a 2 a 



. d_ fdF\ 
Jt urn) 



\-a i R- 



d f dF 



R + 



6k\ dF 



dt \d(DR) 

A remark is necessary at this point. We can also take into account 

Ai =a 3 



a 2 / dR 



OR 



dF dF 
dR + d(VR) 



(2.6) 
(2.7) 



dF 
d(nR) 
(2.8) 



(2.9) 



as a Lagrange multiplier |I3]. The Lagrangian which comes out differs from ( |2.8|) just 
for a term vanishing on the constraint, being 



£ = C-a 3 {R + Q 



- i-) + 

a) 



□ - 



dF 



d(DR) 



(2.10) 



^From this point of view, considering the point-like Lagrangian £ or £ is completely 
equivalent (this remark is obvious dealing with the equations of motion). 
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Let us now derive the Euler-Lagrangian equations from ( |2.8| ). They can be also 
deduced from the Einstein equations ( |2.2|) . For the sake of clarity, let us derive them 
step by step. The equation 

ddC dC 



gives 



+ 2 



dF dF 
dR + nR d(^R) 
d 2 F 



dt da da 



+ 2 3H 2 + 2H + — 



dF 
dR + 



DR-HR 



dR 2 



+ R(aR) 



d 2 F 
d(DRf 



+ 



■ d 3 F ■ 2 d 3 F ■ 

2R 2 ^ + 2{UR) —^— + AR{UR) 



2U 2 R-2H{UR) + R 2 
d 3 F 



d 2 F 



dRd(OR) 



dR 3 
The equation 

gives 
Finally, 



dRd(nR) 2 



dR 2 d(nR) 



0. 



d_dC _ dC 
dtdR ~ dR 



□ - 



0. 



d(nR) 
d dC dC 



dtd{nR) dnR' 
coincides with the Lagrange constraints 

OR = R + 3HR, 
R=-6 (h + 2H 2 + K 



Here H = a/a is the Hubble parameter. The condition on the energy 

C = 0, 



._, . dC -dC . • . dC 

E c = a— + R^ + (DR) ^- 

da dR d{DR) 



which is the (0, 0)-Einstein equation, gives 



dR) 



d fdF\ r 
dt [dRJ + 6 



where 



' 6k\dF dF -d( dF 

R+ ^)dR + DR d(^Rj ~ F - R Jt {d(BR) 



can be interpreted as a sort of effective density (see also |33fl ) . 



+ 

(2.12) 

(2.13) 
(2.14) 

(2.15) 

(2.16) 
(2.17) 

(2.18) 



(2.19) 



(2.20) 
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3 Noether Symmetry Approach 



A Noether symmetry for the Lagrangian ( |2.8|) exists if the condition (|1.5|) holds. It is 

nothing else but the contraction of a Noether vector X, defined on the tangent space 

TQ = {qi,<ii} of the Lagrangian C = £(g^ with the Cartan one-form, generically 
defined as 

ec = ^dq\ (3.1) 
dqt 

Condition (|1.5| ) gives 

iA = S , (3.2) 



where ix is the inner derivative and So is the conserved quantity |21], |22], |23], |35], |36| . In 
other words, the existence of the symmetry is connected to the existence of a vector field 

where at least one of the components a l (q) have to be different from zero. In our case, 
the tangent space is 

TQ = {a, a, R, R, OR, (OR)} , (3.4) 
and the generator of symmetry is 

d d d d-d d 

X = a— + (3— + 7 -r- — + a— + ff— r- + 7 ^- . (3.5) 

da y dR 'd(OR) da H dR d(DR) 

The functions a, /3, 7 depend on the variables a, R, OR. A Noether symmetry exists if at 
least one of them is different from zero. Their analytic forms can be found by expliciting 

( |1 . 5|) , which corresponds to a set of 1 H partial differential equations given by 



2 

2 

equating to zero the terms in a 2 , R 2 , (OR) , aR and so on. In our specific case, n = 3 
is the dimension of the configuration space. We get a system of seven partial differential 
equations 

dF ( da\ d 2 F / dp\ d 2 F ( fry\ 

d 2 Fda d 2 F ( d(3\ d 3 F d 3 F <9 7 d 2 F 

~ 6 dR 2 dR + dRd(nR) [ 3a + 2a dR ) +Pa dR 2 d{uR) +7 ° ' dRd{uRf +a ~dR <9(ni?) 2 ~ ° ' 

(3.7) 

d 2 F da _ d 2 F dp 
dRd(nR) d(DR) ~ a d{DRy d(DR) ~ ' 

dFda d 2 F ( da d(3\ d 2 F ( «9 7 dp\ 

12 mm + 'dR 2 [ 2a + a d^ + a m + a dRWR) [dR ~ 2a d^r 
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2 <9 7 d 2 F n d 3 F d 3 F 
a -z— - „ + 6pa^-^ + 67a- 



dF da 

12 — TTT^rr -f 



9a 0(Di2) 2 
<9 2 F 



dR 2 d(nR) 



(3.9) 



dRd(\JR) 8Rd(nR) 



da d'j 
12a + 6a- + 6a— 



+ 6a 



— a 



_d^F_dp_ 
d(OR) 2 0a 



+ 67a 



<9 3 F 



d 2 F / <9 7 <9/3\ 

3a + a— -, — — + a—— + 



dRd(DR) 2 
d 2 F i 



Q(3a 



d 3 F 



dR 2 d(nR) 



d 2 F dp 
OR 2 d(DR) 



0, 



(3.10) 



d(aRf 



d(DR) OR J dRd(aR) 
d 3 F 



+ Pa 



(dF\ 
F - R— J - paR 



d 2 F 



dRd{URf 
dF 



+ 7a 



( dp da" 
V a d(DR) dR / 

d 3 F 



d 2 F da 
- 6 - — — - + 



OR 2 dUR 



3aOR- 



d(OR)3 
8 2 F 



dR 2 d(nR) 
6k " 



7aDi? 



d 2 F 



dF a d 2 F 

a^z + Pfl77777 + 7 a ~ 



dR 



dR 2 



d{URf dRd{UR) 
d 2 F Q 



(3.11) 



(paR + 7 R)- 



(3.12) 



dRd(aR) 

The system is overdetermined and, if solvable, enables one to assign a, P, 7, and F(R, OR). 
In such a case, we can always transform the Lagrangian ( |2.8|) so that 

C(a, a, R, R, DR, (DR)) -> £(u, u, 117, it>, i) , (3.13) 

where z is a cyclic variable and the dynamics is simplified. This change of variables can 
be easily obtained by the conditions 



dz dz dz 

lxdZ = a d-a +(3 dR +1 dnR = 1 > 
dw dw dw 

lxdw = a ^ +p dR +1 dm = 

du du du 

lxdu = a d^ + p dR + 7 mR = ^ 



(3.14) 
(3.15) 
(3.16) 



which strictly depend on the form of a, f3, 7. Once we solve the dynamics in the system 
{z, w,u}, by the inverse transformation 



{z(t),w(t),u(t)} -> {a(t), R(t), OR(t)} , 



(3.17) 



we recover the dynamics in our primitive physical variables. Here t is the cosmic time. 

However, we have to stress that we are considering a constrained dynamics since the 
variable a, R, OR are related each other. In next sections, we show that the existence 
of the Noether symmetry, i.e. solving the system (|3.6| )- (|3.12|) , gives models of physical 
interest. 
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4 Fourth-Order Gravity 



If the Lagrangian ( |2.8|) does not depend on Di2, we are dealing with fourth-order equa- 
tions of motion. Furthermore, if F(R) = R + 2A, the standard second-order gravity is 
recovered, being A the cosmological constant. 

The configuration space is two-dimensional so that system (|3.6|) - (|3.12|) reduces to 



dF ( da 
dR \ da 



+ a 



dR 2 



P + a 



d£ 
da 



U 2 F\ da 

[dR 2 ) dR ~ ' 

'dF\ da d 2 F 
+ 



dR) dR dR 2 
dF\ 



da dp' 
2a + <X Ta +<X dR. 



0, 



d 3 F 
+ = . 



d 2 F 6k 



dF 



3a F - R— - af3R—— a— + ap 



dR) 



dR 2 



dR 



dR 3 
dR 2 







(4.1) 
(4.2) 
(4.3) 
(4.4) 



where we need only four equations. Immediately we see that a Noether symmetry exists 
for F(R) = R + 2A. Discarding this trivial case, we get the solution 



a 



a 



R 



2P -, F{R)=F R 3 ' 2 



(4.5) 



for any value of the spatial curvature constant k — 0, ±1; Pq and F are constants. This 
case is interesting in conformal transformations from Jordan frame to Einstein frame 
P?| , |3"8[ since it is possible to give explicit form of scalar field potential. In fact, if 



9af3 



'dF\ 

dR) 



In 



'dF\ 
dR) 



(4.6) 



we have the conformal equivalence of the Lagrangians 
C = ^F Q R 3 / 2 < — >£ = 



R 1 




(4.7) 



in our physical units. This is the so-called Liouville field theory and it is one of the few 
cases where a fourth-order Lagrangian can be expressed, in the Einstein frame, in terms 
of elementary functions under a conformal transformation. Using Eqs. ([3.14J) and Q3.15 ), 
we get the new variables 

2 a 2 
w = a R, z = — , 

from which Lagrangian (|2.8|) (without terms in OR) becomes 



(4.8) 



9P ZW ^ 1 /-J 

L = p= — yky/w vtc- 

'w 2 



(4.9) 



A further change of variable y = \Jw gives 



3 

C = 9(3 zy-9ky- 1 L, (4.10) 



where z is cyclic. The dynamics is given by the equations 



11 = 0, (4-11) 

9(3 z + 9k + -y 2 = 0, (4.12) 

3 

9f3 zy + 9ky + ^ = 0, (4.13) 



whose general solution is 



y(t) = Sot + yo, (4.14) 
z (t) = c 4 t 4 + c 3 t 3 + c 2 t 2 + at + c . (4.15) 

Going back to physical variables, we get the cosmological solution 

a(t) = \fWo [cit A + c 3 t 3 + c 2 t 2 + at + c ] 1/2 , (4.16) 

R(t) = (ggj + Vo) 2 u i 7 ) 

1 J 2/5 [c 4 t 4 + c 3 t 3 + c 2 t 2 + at + c ] ' 1 ' J 

The constants q are combinations of the initial conditions. Their values determine the 
type of cosmological evolution. For example, C4 7^ gives a power law inflation while, if 
the regime is dominated by the linear term in a, we get a radiation-dominated stage. 
The system fl4.1| )-( fOD admits other solutions. Another interesting one is given by 

a = Q, P=—, F(R) = F\R + F 2 R 2 , (4.18) 
a 

if the condition 

R=~ (4.19) 

a 1 

is satisfied. Immediately we get the cosmological solution 

a(t) = a t 1/2 . (4.20) 

is 



A similar situation works any time that d 2 F/dR 2 = F" 7^ 0. The system (O 
solved by 

« = 0, /?=Jj. (4-21) 
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Condition ( f4.19| ) has to be satisfied and the radiative cosmological solution ( |4.2CI| ) is 
recovered. Particularly interesting, from the point of view of one-loop corrections of 
quantum gravity, are polynomial Lagrangian of the form 

N 

F(R) = FjR j , N > , (4.22) 

3=0 

where Fj are constant coefficients whose physical dimension depends on j. Cosmological 
models coming from such theories have been widely studied (see e.g. []39, £E|]). However, 
standard Einstein coupling is recovered if F\ = —1/2. 
Results are summarized in Table I. 



Sixth and Eighth— Order Gravity 



Considering OR as a degree of freedom means to take into account the whole system 
( |3.6|) -( j3TT2| ) in order to find some Noether symmetry. If F(R, OR) depends only linearly 
on OR, we have a sixth-order theory, otherwise we are dealing with eighth-order theories. 
A simple sixth-order solution of Noether system 



a: 



fa 



)- (|3.12|) is recovered if 
/?7any, F(R, OR) = F\R + F 2 UR , k = . 



(5.1) 



Due to the above considerations on pure divergence JIT , 12| , this theory reduces to the 
Einstein one where standard cosmological solutions are recovered. 
For powers of OR, we have Noether symmetries given by 



a 



0, P = P , 7 = 0, F(R,OR) = F 1 R + F 2 (OR) n , n > 2. 



However, the theory can assume different forms integrating by parts [0], [12|, |13[ 
equations of motion are 



4FW + QFiH 2 + 2F X — - F 2 (l - n)(DR) 

a 1 

-F 2 n(n - l)a 3 (aR) n - 2 (DR) = S , 
OR — R — 3HR = , 
k 



IV 

a . 







6F X 



F 2 (l - n){UR) n + i?S = . 



(5.2) 
The 

(5.3) 

(5.4) 
(5.5) 

(5.6) 



So is the Noether constant and R is the cyclic variable. The standard Newton coupling 
is recovered, as usual, for F\ = —1/2. A solution for this system is 



a(t) = a$t , 

for k — — 1, So = 0, and for arbitrary n, F 2 . Another solution is 

a(t) = a t 1/2 , 



(5.7) 
(5.8) 
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for k = 0, Eq = 0, F\ = and for arbitrary n, F2. Finally we get 



a(t) = a exp(k t) , 



(5.9) 



for k = 0, S = 0, Fx = 0. 

The radiative cosmological solution ( 4.20| ) with the condition ( 4.19|) is recovered for 
the cases 



a 



0, (3 



0o 



7 any , F(R, OR) = F X R + F 2 R 2 + F 3 DR . 



and 



a 



0, (3 = 0, 7 



To 
a 



F(R, BR) = FxR + F 2 R 2 + F 3 RDR , 



(5.10) 



(5.11) 



with 70 constant. The second one is of physical interest and the related cosmological 



models have been widely studied |L4], |30 

symmetry is rec 

a = , P = /3o, 7 = fa 



Another Noether symmetry is recovered for 

DR 



R 



F(R, OR) = FxR + F 2 VROR . 



(5.12) 



or simply for 

F(R,nR) = F 2 VRDR. (5.13) 

This case deserves a lot of attention since VrTTr is exactly a part of a3-anomaly |L4j 
which can be recovered from the general analysis of one-loop contributions to gravita- 
tional action Z6]. Such a fact is indicative since it seems that searching for Noether 
symmetries could be a relevant method for constructing an effective action of quantum 
gravity. A straightforward change of variables, given by ( J3.14 ) — (3.16), is 



R. 



u 



'DR 
~R~ 



w 



(5.14) 



Choosing the standard Einstein coupling Fx 



-1/2 the Lagrangian (|2.8| ) becomes 



3[ww 2 — kw] — F 2 



3ww 2 u + 3w 2 wit + 



w 3 zit 



2u 2 



— 3kwu 



(5.15) 



where z (i.e. R) is the cyclic variable. ^From (|5.15|) , it is immediate to derive the 
equations of motion. As above, it is possible to recover the particular solutions 



a(t) = dot , a(t) = dot 1 / 2 , a(t) = aoexp(kot) 



(5.16) 



depending on the set of parameters {So, k, F 2 }. A phase-space view and conformal 
analysis, as that in [3D|, gives the conditions for he onset and duration of inflation which, 
specifically, depends on the sign and the value of F 2 and S . This fact restricts the set of 
initial conditions capable of furnishing satisfactory inflationary cosmology as discussed 
in 



14], p0|. In Table II we report the main results of this section. 
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6 Discussion and Conclusions 



In this paper, we have used the Noether Symmetry Approach in order to study higher- 
order theories of gravity. The existence of a symmetry selects the form of higher-order 
Lagrangian as in nonminimally coupled theories (|1.4J) ; there this technique allows to 
assign the form of the coupling F(<p) and the potential V((p) |2l], [22|. Here we discussed 
theories up to eight-order but it is clear that the method works also for orders beyond. 
The scheme is always the same: i) if a symmetry exists, ii) the form of the effective 
Lagrangian is assigned, in) a suitable change of variables allows to write the dynamics 
so that a cyclic coordinate appears. The solution of the cosmological problem results 
simplified since two first integrals of motion are present (the energy and the symmetry). 
However the scheme works if a suitable minisuperspace has been a priori defined. Here 
we used FRW minisuperspaces. 

Some final remarks are necessary at this point. First of all, by this technique, one is 
capable of selecting higher-order terms as i? 3//2 or y/ROR of physical interest since they 
can be connected to the one-loop or trace anomaly corrections of the effective action 
of quantum gravity. It is worthwhile to stress that such terms are not perturbatively 
introduced but emerge by the request of symmetry. Furthermore, the system of partial 



differential equations (|3.6|) - (|3.12|) can have several solutions and their finding out is just 
a question of mathematical ability. In this paper, we have not made an exhaustive list 
of the possible Noether symmetries in higher-order theories, but we have only presented 
some examples in fourth-, sixth-, and eighth-order models. 

A further point which has to be stressed is that obtaining some higher-derivative 
terms in the effective Lagrangian of gravity does not, automatically, makes the theory 
renormalizable. In other words, functions of R and OR alone do not give renormalizable 
Lagrangians. For renormalization in fourth-order gravity, see e.g. f7|. Our point of view 
is that the existence of a Noether symmetry for the Lagrangian seems to be connected to 
the existence of corrective terms which one needs for renormalization and the approach 
presented in this paper seems a scheme which could be generalized. 

Finally, the existence of a Noether symmetry makes the analysis of a given cosmology 
more tractable; however the existence of such a symmetry is not, in itself, a sufficient 
motivation for preferring a particular theory. The situation, as also shown elsewhere 
22) , p3[ , 24j| , becomes interesting if the selected theory is physically relevant per se and 



the Noether approach selects just it. 

A further step which the authors are going to achieve is to apply the Noether Sym- 
metry Approach to the full field theory without reducing to particular minisuperspaces 
as done here. 
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Table I- Symmetries in Fourth-Order Models 



F(R) 


a 


p 


R + 2A 





P(a,R) 


F RV 2 


Pod' 1 


-2(3 a- 2 R 


F R + FxR 1 







TUF W 





PoiaF")- 1 



Table II- Symmetries in Higher than Fourth-Order Models 



F(R, OR) 


a 


p 


7 


FiR + F 2 (OR) n ( n ^ 1) 





Po 





FiR + F 2 R 2 + F 3 OR 





P a~ l 


7(a,i2, OR) 


F R + F X R 2 + F 2 ROR 








7 a _1 


F 2 ^R{OR) 





Po 


PoR^OR 


F 1 R + F 2y jR(DR) 





Po 


PoR^OR 
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